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Abstract 

We consider a regular random network where each node has exactly 
three neighbours. Ising spins at the network nodes interact antiferromag- 
netically along the links. The clustering coefficient C is tuned from zero 
to 1/3 by adding new links. At the same time, the density of geometri- 
cally frustrated links increases. We calculate the magnetic specific heat, 
the spin susceptibility and the Edwards- Anderson order parameter q by 
means of the heat-bath Monte Carlo simulations. The aim is the tran- 
sition temperature Tx dependence on the clustering coefficient C. The 
results are compared with the predictions of the Bethe approximation. 
At C = 0, the network is bipartite and the low temperature phase is an- 
tiferromagnetic. When C increases, the critical temperature falls down 
towards the values which are close to the theoretical predictions for the 
spin-glass phase. 
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1 Introduction 

Statistical mechanics of random networks gained recently many applications in 
interdisciplinary sciences. The list of references is already very rich; in almost 
each year new monographies appear [U [U [31 |H [S] . In theory, the promising 
challenge is to investigate collective phenomena in networks . More than often 
it is useful to decorate nodes i with additional variables, as Ising spins Si = ±1. 
Properties of such networks are of interest for theory of computations, as infer- 
ence problems, but also for the theory of disordered magnetic systems. Here we 
are interested in the spin-glass phase of the small-world random networks. 

For trees, where closed loops are absent, there is a consistent theory of Ising 
magnetism, i.e. the Bethe theory [8l[9j[T0]. Once the loops appear in the sys- 
tem, the Bethe theory becomes an approximation. In our system, small loops 
are introduced when we enhance the clustering coefficient [TTl [T^]. Still, the 
results of the Bethe theory are useful as a point of reference. As it was demon- 
strated in [I^, the accordance of the numerical results with this theory was 
deteriorated when C increases. On the other hand, there is at least one periodic 
two-dimensional Ising lattice when all bonds are antiferromagnetic, the density 
of frustration is high and the transition temperature is positive; this is one of the 
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Figure 1: The magnetic susceptibility x calculated from the variance of the mag- 
netization against temperature T, for different values of the clustering coefficient 
C=0.0, 0.08, 0.f4, O.fS, 0.25, 0.29 and 0.33. With increasing C, the maximum of 
X shifts to low temperatures till C=0.f8, then remains approximately constant 

Archimedean lattices [H]. The ground state of this system is highly frustrated, 
but at least some of the energy barriers between different states remain finite. 
Although at zero field the net magnetization is zero, there is no disorder in 
the system; therefore it is hard to speak about the spin-glass phase. Then, the 
(3, 12^) Archimedean lattice can be compared to our system in the case where 
C is maximal. 

Recently we investigated the transition from the paramagnetic to the spin- 
glass phase in the random Erdos-Renyi network with enhanced clustering co- 
efficient C [T^. The enhancement is introduced by adding new links between 
neighbours of the same sites [TT]. In this way, C varied from almost zero to 
about 0.3. The main result was the transition temperature Tsa dependence on 
the clustering coefhcient C. However, in the numerical plots the transition was 
partially hidden by the contribution of spins which could flip with zero energy 
cost. That is why here we investigate the same transition in the regular random 
network, where the degree of each node is an odd number - here it is equal 
to three. In such systems no spin can flip at zero energy cost in zero external 
magnetic field. 

In the next section we show the method and the numerical results. In Sec- 
tion III the transition temperature is compared with the predictions of the 
Bethe theory. Short discussion closes the text. 
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2 Method and results 



Initial form of the constructed network is prepared as an even number iVo of un- 
linked nodes. The system is divided into two equal parts. Three 'dangling' links 
are assigned to each node. Then, the links of one part of nodes are randomly 
joint to links of the second part. In this way, initially the system is bipartite, 
with the clustering coefRcient C equal to zero provided that the number of nodes 
is large enough. Next step is performed for each node with the probability p: 
the node is substituted by a group of three nodes, mutually linked. After this 
substitution, three previous neighbours of the node arc linked to the nodes of 
the group. In this way, the degree of each node remains exactly three. As the 
Ising interaction is exclusively antiferromagnetic, the frustration is purely geo- 
metrical. When C = 0, there is no frustration. The density of frustrated bonds 
increases with C . li p — \, the clustering is maximal and for each node, two 
out of its three neigbours are linked to each other; then C — 1/3. The final 
number of nodes is iV = iVo(l + 2p); for each p the amount iVo is chosen as to 
get approximately the same N . 

The heat-bath Monte Carlo algorithm is applied to investigate the magnetic 
properties: the susceptibility x, the specific heat Cy and the Edwards- Anderson 
order parameter q. x can be calculated from the field derivative of the magne- 
tization m{h) or from the variance of the spectrum of m(0). Similarly can 
be obtained from the temperature derivative of energy or from the variance of 
the energy fluctuations. The parameter q is calculated from the formula [14j 



The antiferromagnetic exchange integral J is set to -1. As a rule, the calcula- 
tions are performed for N equal at least 9 x 10^ spins. The time of calculation - 
after a necessary transient - was usually 10^ steps, where one step is equivalent 
to probing N spins in random order. The results are limited to temperature 
T > 0.5, where the transient time for relaxation of energy was smaller than 10'^ 
timesteps. 

In Fig. 1 we show the magnetic susceptibilities xi^) for different values of 
the clustering coefficient C. The temperature is expressed in energy units — J. 
The finite size effect is checked for C = and C = 0.33; the differences between 

= 9 X 10^ and A = 9 x lO'^ are invisible. As a rule, the plots obtained from 
the field derivative of the magnetization coincide with those from the variance 
of the magnetization at zero field above the transition temperature. Below r^,, 
the curves split for intermediate values of C (between 0.1 and 0.29). The plots 
for the magnetic specific heat obtained from the variance of energy (Figs. 2 and 
3) and from its thermal derivative coincide in most cases, even below Tx and for 
intermediate values of C . However, the data from the maxima of the specific 
heat are known to provide an evaluation of the upper limit of the transition 
temperature rather than itself. The results on obtained by different 
methods are compared in Fig. 5. As we see, there is a systematic split between 
Tx obtained from x ^^id obtained from C„, but the character of the curve 
Tx{C) is preserved. In Fig. 4 we show the data on obtained from the thermal 
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Figure 2: The magnetic specific heat against temperature T for different val- 
ues of the clustering coefficient C—Q.O (squares), 0.08 (stars) and 0.14 (circles). 
With increasing C, the maximum of Cy shifts monotonously to low tempera- 
tures. 



dependence of the Edwards- Anderson order parameter q. These data also show 
a decrease of with C. However, this decrease persists till T = 0.5, where our 
numerical results are less reliable. 



3 Discussion 

For C = the network is bipartite. Although the interaction is antiferromag- 
netic, the system is equivalent to a ferromagnet. To see this, we have to flip 
spins at half of nodes and to invert simultaneously the sign of all exchange in- 
tegrals J from -1 to The formula for the Curie temperature (J = -1-1) for 
the regular Bethe lattice with the coordination number A; is [1] 

Tc = j^ (2) 

In this case this formula can be applied equally well to the Neel temperature 
(J — —1). In both cases the number of loops in the system is exactly zero. 

For C > the system contains loops and it is no more equivalent to a 
ferromagnet. On the other hand, the network is no more a regular Bethe lattice. 
The number of second neighbours decreases gradually with p from 6 when p = 
to 4 when p — 1. The latter case is, however, in some sense more regular 
than the one with smaller p, because all nodes have again the same number 
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Figure 3: The magnetic specific heat against temperature T for different 
values of the clustering coefficient C= 0.18, 0.25, 0.29 and 0.33. With increasing 
C, the height of the maximum of decreases. 



of the second neighbours. Then the Neel temperature (J = —1) between the 
paramagnetic and the antiferromagnetic phase is to be found from 

B-l 

and the transition temperature between the paramagnetic and the spin-glass 
phase - from 

where B is the average branching parameter, i.e. B = 22/^1, where zi (22) is 
the average number of first (second) neighbours ^ . 

Although the data from the specific heat fit better to the Bethe theory, it 
is known that the position of the maximum of is rather the upper bound of 
the transition temperature than this temperature itself. 

The numerical results indicate, that for C = the low-temperature phase is 
antiferromagnetic. They also suggest that for C > 0.1 and T < the system 
becomes a spin glass. However, the identity of this phase remains not clear. For 
C = 1/3 the system can be compared to the (3, 12^) Archimedean lattice [13] 
with random rewiring. We can expect that the ground state is degenerated in 
the same way as in this Archimedean lattice. However, the condition is that the 
rewiring does not destroy the small triangles introduced when C is enhanced. 
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Figure 4: The Edwards- Anderson order parameter q against temperature T for 
different values of the clustering coefficient C=0.0, 0.08, 0.14, 0.18, 0.25. With 
increasing C, the temperature Tx where q vanishes shifts monotonously to low 
temperatures. 



In this way the frustration is not altered by rewiring: the frustration remains 
local. Then, there is frustration and disorder, but the disorder does not influ- 
ence the frustration. On the other hand, the energy barriers between the local 
energy minima remain finite, because the system can move from one minimum 
to another by flipping a finite amount of spins. Although the system behaves 
as a spin glass, some ingredients of this mysterious phase are missing. 

To summarize, our intention was to compare the numerical results on the 
transition temperature with the results of the Bethe approximation. In Ref. 
[Hj , distinct departures have been found between the Bethe theory and the nu- 
merical experiment. Here the departure is smaller and quantitative rather than 
qualitative. As a rule, the transition temperature is overestimated by the Bethe 
approximation. Most important difference between the system considered in 
this text and the system discussed in Ref. [12] is the degree distribution. This 
suggests, that the qualitative departure of the results of Ref. from those of 
the Bethe approximation is due to the variance of the degree distribution. 
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Figure 5: The critical temperature against the clustering coefficient C, cal- 
culated numerically from the inflection point of the magnetic susceptibility 
X (stars), the maximum of the specific heat Cv (squares) and the Edwards- 
Anderson order parameter q (circles). Two lines mark theoretical values of 
for the para-antiferromagnetic phase transition (Eq. 3, upper line) and the 
paramagnetic- spin glass transition (Eq. 4, lower line). 



7 



References 

[1] A.-L. Barabasi, Linked: How Everything Is Connected to Everything Else 
and What It Means for Business, Science, and Everyday Life, Plume Books, 
New York 2003. 

[2] S. N. Dorogovtscv and J. F. F. Mendes, Evolution of Networks: From 
Biological Nets to the Internet and WWW, Oxford UP, Oxford 2003. 

[3] Handbook of Graphs and Networks: From the Genome to the Internet, Eds. 
S. Bornholdt and H. G. Schuster, Wiley- VCH, Berlin 2003. 

[4] R. Pastor-Satorras and A. Vespignani, Evolution and Structure of the In- 
ternet: A Statistical Physics Approach, Cambridge UP, Cambridge 2004. 

[5] R. Durrett, Random Graph Dynamics, Cambridge UP, Cambridge 2006. 

[6] G. Caldarelli, Scale-Free Networks, Oxford UP, Oxford 2007. 

[7] S. N. Dorogovtsev and A. V. Goltsev, Critical phenomena in complex net- 
works. Rev. Mod. Phys. 80 (2008) 1275. 

[8] S. Katsura and M. Takizawa, Bethe lattice and Bethe approximation. Prog. 
Theor. Phys. 51 (1974) 82. 

[9] R. J. Baxter, Exactly Solved Models in Statistical Mechanics, Academic 
Press, London 1982. 

[10] J. M. Mooij and H. .J. Kappen, On the properties of the Bethe approximation 
and loopy belief propagation on binary networks, J. Stat. Mechanics: Theory 
and Experiment (2005) P11012. 

[11] P. Holme and B. J. Kim, Growing scale-free networks with tunable cluster- 
ing, Phys. Rev. E 65 (2002) 026107. 

[12] A. Manka, K. Malarz and K. Kulakowski, Clusterization, frustration and 
collectivity in random networks, Int. J. Mod. Phys. C 18 (2007) 1765. 

[13] M. J. Krawczyk, K. Malarz, B. Kawecka-Magicra, A. Z. Maksymowicz and 
K. Kulakowski, Spin-glass properties of an Ising antiferromagnet on the 
Archimedean (3, 12^) lattice, Phys. Rev. B 72 (2005) 24445. 

[14] K. Binder and A. P. Young, Spin glasses: Experimental facts, theoretical 
concepts, and open questions. Rev. Mod. Phys. 58 (1986) 801. 



8 



